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Abstract 

A model-free theoretical framework for a phenomenological description of spin- 
lattice relaxation by anomalous translational diffusion in inhomogeneous systems 
based on the fractional diffusion equation is developed. The dependence of the spin- 
lattice relaxation time on the size of the pores in porous glass Vycor is experimen- 
tally obtained and found to agree well with our theoretical predictions. We obtain 
nonmonotonic behavior of the translational spin-lattice relaxation rate constant (it 
passes through a maximum ) with the variation of the parameter referring to the 
extent of inhomogeneity of the system. 
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1 Introduction 



Anomalous diffusion processes in randomly disordered media are of consider- 
able interest at present (see [1] and refs. therein). The list of examples includes, 
glasses and supercooled liquids [2], porous [3], [4], [5], percolative [6], [7], [8], 
polymeric [9] and diffusive [10] systems, etc. One of the powerful methods for 
investigation of such processes is that of nuclear magnetic resonance (NMR) 
diffusometry [3] . The latter was successfully applied to the systems listed above 
[3-9]. The theory of NMR diffusometry in disordered media [4], [8] is devel- 
oped within the framework of a modern approach to anomalous diffusion [1] 
invoking the so called fractional calculus [11]. 

The aim of the present paper is to report the experimental results on spin- 
lattice relaxation time in porous glass Vycor with molecules of hexane as a 
diffusion tracer and to compare them with phenomenological generalization 
of the theory of spin-lattice relaxation in homogeneous media for the case of 
inhomogeneous one. 

A modern description of anomalous diffusion processes in disordered (e.g., 
porous) media makes use of the so called fractional diffusion equation (FDE). 
The idea of this approach goes back to many pioneers whose achievements are 
honored in the review article [1]. The reason for introducing this equation is 
as follows. As is well known the mean squared displacement of a free particle 
in a homogeneous media increases linearly with time < x^{t) > ~ Cit where 
Ci is the diffusion coefficient with the dimension cw?/s. This conventional 
Einstein relationship of the classical theory results from the ordinary diffusion 
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equation for the probability density function to find the particle at position 
X at time t and is a direct consequence of the Pick's second law. In the case 
of an inhomogeneous media a disorder (like e.g., obstacles caused by walls of 
pores in a porous materials) leads to a slower increase of the mean squared 
displacement with time 



where Ca is the diffusion coefficient with the dimension err? js"" and < a < 1 
is a phenomenological parameter characterizing the extent of inhomogeneity. 
Such behavior is called sub-diffusion and it originates in any fractal media due 
to the presence of dead ends on current ways. It should be stressed that the 
departure from the classical theory is of fundamental character and the latter 
relationship can not be reduced to the former one by simple scaling of time and 
diffusion coefficient. The most straightforward way to obtain mathematically 
sub-diffusion is to generalize the ordinary diffusion equation by replacing the 
ordinary derivative in time by the fractional one of the order a. The resulting 
FDE is a convenient mathematical tool although a phenomenological one since 
there is no lucid and commonly accepted physical meaning of a fractional 
derivative at present. Some theoretical motivation for introducing the FDE 
arises from its intrinsic relationship with the so called continuous time random 
walk theory. In fact the FDE can be derived from a generahzed Langevin 
equation with the memory kernel accounting for a power-law waiting time 
statistics of the trapping events [1] . 

Much work is being carried out at present to verify experimentally the predic- 
tions of the approach based on the FDE [4], [8], [10]. As is emphasized in the 

papers [4], [8] NMR microscopy enabled the authors to verify the anomalous 



< X^{t) > ~ Cat' 
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solutions of the FDE for the first time. The results of these work are com- 
patible with other approaches to NMR relaxation in porous media [28], [29]. 
However the fact that the authors of [4], [8] resort to a one-dimensional frac- 
tional counterpart of the ordinary diffusion equation makes their theoretical 
calculations essentially of qualitative character and leaves room for further 
advances. The present paper is a development along this line. We resort to a 
three-dimensional FDE and obtain stringent quantitative relationship of the 
spin-lattice relaxation rate constant with the parameters (a and C^) of the 
FDE. Our results provide new and useful analytical tool to interpret the man- 
ifestation of anomalous translational diffusion in NMR relaxation which is 
substantiated by the obtained experimental data. 

The matter of NMR relaxation at diffusion through a pore is of particular 
significance for biological systems. According to a widely accepted point of 
view one of the ways of molecule penetration through a biological membrane 
is associated with the short-lived and long-lived pores with a diameter of 
several Engstroms [12]. Such pores are identified with proteins - aquaporins 
revealed in cellular membranes of various organisms [13], [14]. In that case the 
water molecules can move only one by one according to the so called single 
file diffusion. Such pores seem to be mostly suitable for providing selection 
and thin regulation of the transfer but the question about their reality and 
peculiarities of their functioning still remains a matter of discussion. 

At present the dynamics of molecules of liquids in porous systems attracts con- 
siderable attention with an accent on the questions of multi-dimensionality of 
the diffusion and the existence of bound motion [15], [16]. The diffusion of 
water in pores with a diameter comparable with that of water molecules is 
studied by molecular dynamics [17-20]. In [18] it is shown that fluid mechan- 
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ics (classical hydrodynamics) can be qualitatively extrapolated for the one by 
one motion of the water molecules in pores with the diameter of 3 Engstroms. 
Water in a narrow pore retains the properties of liquid state with the diffu- 
sion coefficient of approximately 70% of that in the bulk phase. In [19], [20] it 
is shown that water in narrow pores has the properties of a structured liquid 
with angular and radial distributions and with the diffusion coefficient compa- 
rable with that of the bulk phase. According to the data of neutron diffraction 
spectroscopy the characteristic time of passing of the water molecule through 
the membrane is of order of 100 mks [21] and thus finds itself within the time 
scale of the processes amenable to NMR spectroscopy. It is reasonable to ex- 
pect that the constraint of the water molecules' motion and their interaction 
with the walls of the pores leads to peculiarities of magnetic relaxation of 
water. Our final interest in transport processes in biological membranes mo- 
tivates the present attempt to consider magnetic relaxation by translational 
diffusion in the pores of much simpler object namely a porous glass. 

Concluding the Introduction we would like to emphasize the following fact. 
The FDE is much more adequate than the ordinary diffusion equation for the 
description of processes in inhomogeneous media but one has to pay a hard 
cost for gained facilities. Namely the diffusion coefficient in the FDE acquires 
a functional dependence on the phenomenological parameter a referring to the 
extent of inhomogeneity. This dependence is not known in general case and 
at present one can only guess some trial possibilities for particular systems 
of interest. Besides this diffusion coefficient acquires the unusual dimension 
cm^/s" and thus has minimal crossing (only at a = 1) with the ordinary 
diffusion coefficient with the dimension err? js which can be obtained by direct 
measurements employing the existing scheme of Pulsed-Field Gradient NMR. 
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2 Results 



2.1 Theory 

The theory of spin-lattice relaxation by dipole-dipole interaction initiated by 
Bloembcrgcn, Purccll and Pound and developed by Torrey and Abragam is 
presented in [22]. The rate constant of spin- lattice relaxation 1/Ti is decom- 
posed into translational and rotational parts 

l/Tl = {1/T,)trans + {l/T,)rot (2) 



To calculate this value theoretically one needs to know the spectral density 
of the correlation function for spherical harmonics. The relationship of the 
value of the spectral density at a Larmor frequency ujl with the value of the 
contribution to the spin-lattice relaxation rate constant due to translational 
diffusion is of the form 

{m)trans = ^T^^'^U + 1) ^^(1) (^^^ ^ j(^)(2u;^)} (3) 



where 7 is the gyromagnetic ratio of the nucleus, / is their spin and h is the 
Planck constant. The spectral densities are proportional to the function J{uj) 
of the spectral density: J^^\uj) = ff and J^'^\uj) = ^J{(^)- Thus one 
actually needs to calculate the only function J{uj). The consideration of the 
spin-lattice relaxation in homogeneous media due to diffusion motion of the 
particles presented in [22] is based on the ordinary diffusion equation for the 
probability density. 
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How is the value of 1/Ti altered when we pass from a homogeneous system to a 
disordered (e.g., porous) one ? One can expect that {l/Ti)trans is very sensitive 
to the extent of porosity in the whole range while {l/Ti)rot becomes to feel 
the porosity only when the size of the pores becomes commensurate with that 
of the molecules, i.e. in a wide range of porosity it remains constant. That is 
why namely the dependence of the value {1/Ti) trans on the extent of porosity 
is of interest for theoretical consideration in the present paper. We develop the 
latter within the framework of a modern approach to anomalous translational 
diffusion in disordered systems based on the FDE. The reason to resort to 
such generalization of the ordinary diffusion equation for the description of 
anomalous diffusion processes is stated in the Introduction. 

The generalization of the ordinary diffusion equation within the fractional 
calculus was suggested by Schneider and Wyss [23] 

^^^cMr^'P)ir,t) (4) 

where P{f, t) is the probability density function to find the particle at position 
f at time t, is the three-dimensional Laplace operator, denotes the 
fractional diffusion constant with the dimension [cm /s"] and D^+" is the 
Riemann-Liouville fractional derivative of order 1 — a and with lower limit 0+ 
which is defined via the following relationship [11] 

^ ^ 

where r(x) is a gamma function. The solution of equation (4) for the case of 
sub-diffusion < a < 1 with the initial condition P(f, 0) = S{f) where 5{x) 
is a Dirac function is obtained in [23] and expressed via the Fox's function 
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[24],[25],[26] 



P{f,t) 



tt20 
^12 



\ 



{l,a) 

(3/2,1), (1,1) 



(6) 



The latter is defined as i^™" 



V 



ib,,B,),...,{bp,Bp) 



2-Ki 



J 



where ryls = -np — , . sf4r — tt^u — s-^- The nomenclature m the Fox s 
function associated with the vertical bar is explained via its explicit definition 
by the contour integral. The requirements to the contour path L are formu- 
lated in [24]. 

Now we follow the algorithm of [22] (in what follows all corresponding results 
from [22] are obtained as a particular case a=\ of the present approach). Under 
f we denote the vector f\ — fi connecting two identical molecules diffusing 
relative to each other rather than the radius- vector of the molecule diffusing 
relative to a fixed point. This leads only to the change of 4Cat° by SCq^" in 
(6). Our aim is to calculate the correlation function 

Tf (^(0),^(0))T-(^(t),(^(t)) 



G{t) = N 



tI 



-P(f -fo,t)d-^rod-^r 



(7) 



where N is the number of spins in 1 cm^, T^{6,(p) is a spherical harmonic 
and * denotes complex conjugate. To be more precise we need the spectral 
density of the correlation function G{t) to calculate the spin-lattice relaxation 
rate constant with the help of (2). At integration in (7) one should take into 
account that r and Tq can not be less than some limit value d - the least 
distance to which the molecules can approach to each other. If the molecules 
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are considered as spheres of the radius a then d — 2a [22] . 



Making the Fourier transforming of the function P(f — fo, t) with the men- 
tioned above change of ACat"" by SCo-t" and hterally repeating the manipula- 
tion of [22] with spherical functions one obtains 

m 7du 



G{t) 











0-20 
^12 



i?2 



V 



(3/2,1), (1,1) 



(8) 



where Ju{x) is a Bessel function of order u and i? = |f — fo|. The integration 
can be fulfilled and yields the expressions for both the correlation function 
and its spectral density via Fox's functions. However a Fox's function is to 
regret not tabulated at present either in Mathematica or Maple or Matlab. 
Thus it is rather difficult to use such formulas for plotting the behavior of the 
correlation function and the spectral density. That is why it is useful to obtain 
another representation of the spectral density which enables one to plot the 
frequency dependence of this function. For this purpose we make the Fourier 
transforming of the Fox function in (8). The latter can be achieved with the 
help of the following trick going back to original investigations of Fox. In [23] 
the Mellin transform of the Fox function is given 



Ml 



I 



tt20 
-^12 



(l,a) 

(3/2,1), (1,1 



1 



a V(8C„)i/2 
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Now we make use of the identity [27] 

Fc {f{x)-uj} = |r(s) cos y M {/(x); 1 - (10) 

where Fc denotes the cosine Fourier transform and M denotes the MeUin 
transform which for our case is given by formula (9). We denote the charac- 
teristic time 



2Cq; 



:iii 



As a result we obtain 

Nt^^"" na f r 

J(uj) — — 1— sin — / du Js (u) 



oo 

r 1 2 

X 





u 



U'^ + {uJTa) " + {oJTa)" COS ■ 



[12) 



The integral in (12) can be easily calculated numerically. We use Mathematica 
for this purpose. Formula (12) is the main theoretical result of the present 
paper. This formula is a direct generalization of the textbook result (see e.g., 
[22] formula VIII. 114) for the homogeneous case a = 1 for the case of arbitrary 
< q; < 1. The spectral density is finite at — > only for a—1. At all a < 1 
it is power- law divergent in this limit. 

2.2 Phenomenology 

The predictions of the present theory are defined by the dependence of the 
diffusion coefficient on the phenomenological model parameter a referring 
to the extent of inhomogeneity. This dependence is of most interest because it 
characterizes the disordered media [1], [8], [10]. Its derivation is a separate and 
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still open problem. It should be a matter of a fundamental theory and is out of 
the scope of the present paper. In our phenomenological approach we restrict 
ourselves by a model example. As such we choose a particular behavior of the 
diffusion coefficient 



where Ci is given by the Stokes formula 

c. = lsL (14) 

bnar] 



Thus the characteristic time is parameterized as 

Ta = rir(a)i/" (15) 



where 



The behavior of the spectral density at several values of the parameter a 
is depicted in Fig.l. Substitution of (12), (13) and (15) into (3) yields the 
normalized spin-lattice relaxation rate constant 

5dCi 



(<^lTi)"~^ r(Q;) sin — j du J^{u) 



oo 

r 1 2 

ux 



[ + iuLnY'' r{ay + 2v? {ult^T T{a) cos ^ 

2^ 

m4 + {2ujLTif'^ T{af + 2^2 {2ujLTif T{a) cos ■ 



(17) 
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Its behavior on the parameter of inhomogeneity a at different values of ojlTi 
is depicted in Fig. 2. At experimentally used Larmor frequency wlTi fa 0.015 
the rate constant of spin-lattice relaxation has a maximum at a ~ 0.35 (with 
WLTa ~ 0.21). 

2.3 Objects of investigation, experimental procedure and results 

Porous glasses of the class "Vycor", with characteristics presented in Table 
1, were used as model samples of porous media. There R is the size of glass 
microparticles, p is the pore radius, Sn and are the pore specific surface 
area and volume, respectively. The primary porosity of these glasses is deter- 
mined by the pore size p, and the secondary one is determined by the glass 
microparticle size R. 

In the experiment, hexane vapor was absorbed at the temperature of 21°C in 
the H-shapcd cell consisting of two connected tubes with the inner diameter 
of 7 mm. One tube was filled with a sample of porous glass to the height of 
12 mm, and into the other tube the liquid hexane was poured. Hexane vapor 
absorbtion took place in the presence of a certain amount of air. Prior to cell 
soldering, the leg with porous glass was heated at the temperature of 140°C 
to remove water vapor. The amount of the absorbed hexane was determined 
from the amplitude of free induction decay after the 90° radio frequency pulse. 
The calibrating of the signal was carried out using the reference sample with 
the known contents of protons. The process of absorption consists of two 
stages. The first faster stage has the character of mono molecular absorption 
and is terminated at the amount of absorbed molecules equal or close to the 
amount necessary to form a mono molecular layer. The second slow stage is 
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associated with the capillary condensation, and lasts from 24 hours ("Vycor 
20") to a month and a half ("Vycor 200"), and is accomplished by fiUing 
the micropores. A good coincidence of calculated and measured values of the 
amount of hexane, necessary to fill completely the primary pores, is observed. 
Spin-lattice relaxation times, Ti, of absorbed hexane were measured using a 
zero-method of the pulse sequence 180'' — 90° on the coherent NMR relaxation 
meter using protons at the frequency of 19.5 MHz. Taking the average size of 
a hexane molecule as a = 5A {d = 2a = lOA) and its diffusion coefficient in 
bulk water as 4.17 * 10~^cm^/s at room temperature we obtain ti 10~^°s 
and consequently ujlTi ~ 0.015. 

The decay of the longitudinal magnetization of the hexane in the region of full 
filling of the pores is one-exponential for all porous glasses. As an example the 
decay of the longitudinal magnetization of the hexane in the "Vycor 160" is 
presented in Fig. 3. The experimentally observed spin-lattice relaxation times 
depending on the inverse radius of the pores are plotted in Fig.4. One can 
see that with the variation of the pore size from 400 A to 55 A the rate 
of longitudinal relaxation increases and for p=20 A it begins to decrease in 
accordance with the prediction of the theory. 

3 Discussion 

Our theoretical approach is based on the FDE like that of the papers [4] , [8] . 
However in these papers a one-dimensional fractional counterpart of the ordi- 
nary diffusion equation is employed to gain qualitative insight. In contrast we 
consider a three-dimensional one that enables us to develop quantitative for- 
malism. In our approach we do not a priory assume any type and character of 
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the motion like e.g., quasi- two-dimensional bulk mediated surface diffusion of 
the absorbed molecules in the pores [4]. Our approach due to its phenomeno- 
logical model-free character is flexible enough to include formally such type of 
motion as a particular case. 

Translational diffusion was previously taken into account within the so called 
reorientation mediated by translational displacements model [3], [4]. It is as- 
sumed there that dipolar coupling responsible for the relaxation mechanism is 
predominantly of an intramolecular nature, i.e. the fluctuations causing relax- 
ation are exclusively due to molecular reorientations [3]. In our approach we 
emphasize another facet of translational diffusion manifesting itself in the in- 
termolecular interactions. We conceive the process as a bulk three-dimensional 
diffusion in space with obstacles imposed by the porous media. In our opinion 
the porosity of any extent considerably manifests itself in the translational 
diffusion. On the other hand it can affect rotational diffusion appreciably only 
when the size of pores becomes commensurate with that of the molecules. 
That is why the value of {Ti)rot is independent on the size of the pores in a 
wide range of the extent of porosity while {Ti)trans is very sensitive to this 
parameter. The latter value characterizes the porosity of the system and is of 
main interest in the present paper. 

The most important feature of the spectral density (12) which distinguishes 
it in the general case from the particular one a = 1 considered in [22] is the 
divergence of the spectral density at low frequencies (see Fig.l). Our phe- 
nomenological model-free approach can be compared with that based on a 
molecular model of the porous system [28]. Fig.l is qualitatively similar to 
Fig. 2 from [28]. However there are significant quantitative distinctions be- 
tween two approaches. The authors of [28] obtain a logarithmic divergence of 
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the spectral density at low frequencies while we obtain a power-law one. Our 
result agrees with the conclusions of the papers [4], [29] where a power- law 
divergence 1/Ti oc at low frequencies is reported. However the authors of 
these papers conclude that the proton relaxation process is mainly due to in- 
tramolecular dipolar interaction while we explore the range of inhomogcneity 
(particularly porosity) of the system where in our opinion intermolecular inter- 
actions of molecules determine the dependence of the relaxation rate constant 
on translational diffusion process. 

It is known that a confinement enhances significantly both spin-lattice and 
spin-spin relaxation rate constants and alters their frequency and tempera- 
ture dependencies (see [28] and refs. to previous works of its authors). In 
these papers a monotonic behavior of the spin- lattice relaxation rate constant 
(1/Ti cx 1/p or 1/Ti cx where p is the average pore size) is reported. 
In contrast we obtain nonmonotonic behavior of the translational spin-lattice 
relaxation time with the variation of the parameter a referring to the extent of 
inhomogeneity ( or that of porosity, namely the size of the pores) with a — 1 
corresponding to the case of homogeneous system. We find that at a given 
temperature there is an optimal porosity for maximal spin-lattice relaxation 
rate constant. Such behavior is determined by contribution from translational 
diffusion. For experimentally used Larmor frequency wi,Ti ~ 0.01 and for a 
particular model dependence of the fractional diffusion coefficient on the pa- 
rameter of inhomogeneity the rate constant of relaxation first increases with 
the decrease of a then passes through a maximum at wiTa ~ 0.1 and finally 
sharply decreases. The reason for such behavior is quite similar to that for or- 
dinary diffusion in a homogeneous system with the variation of temperature. 
In the latter case the characteristic time ti(T) is a function of temperature 
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and the rate constant of spin-lattice relaxation has a maximum at some tem- 
perature where ulTi ~ 1. For inhomogeneous systems the characteristic time 
Ta{T, a) besides temperature is a function of the parameter of inhomogeneity 
a and the variation of the latter (e.g., that of the porosity of the system) at 
constant temperature leads to a maximum of the spin-lattice relaxation rate 
constant at some a. However this maximum in general case exists at coLTa 
different from 1. 

4 Conclusions 

A generalization of the ordinary diffusion equation within the framework of 
the fractional calculus provides natural phenomenological generalization of 
the theory of spin-lattice NMR relaxation in homogeneous systems for the 
case of the inhomogeneous ones. This development is carried out within a 
general modern trend to extend the theory of Gaussian translational diffusion 
in homogeneous systems for the case of anomalous diffusion in disordered 
(e.g., porous) systems by making use of the fractional diffusion equation. The 
translational contribution into spin-lattice relaxation time is found to be highly 
sensitive to the extent of inhomogeneity. It exhibits nonmonotonic behavior 
with the variation of the parameter a referring to the extent of inhomogeneity 
(with a = 1 corresponding to the case of a homogeneous and isotropic system 
and a — corresponding to that of an absolutely inhomogeneous one). The 
rate constant of relaxation passes through a maximum with the decrease of a. 
For experimentally used Larmor frequency wlTi ~ 0.015 and for a particular 
model dependence of the fractional diffusion coefficient on the parameter of 
inhomogeneity the maximum is at a ~ 0.35 (wlTq ~ 0.21). One can conclude 
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that the present work provides a rehable theoretical framework for the analysis 
of anomalous translational diffusion processes by NMR microscopy. 

We obtain the explicit dependence of the rate constant of spin-lattice relax- 
ation by anomalous translational diffusion on the parameter of inhomogeneity 
a taking place in the fractional diffusion equation. Formula (12) is a direct 
generalization of the textbook result (see e.g. formula V11I.114 in [22]) for the 
homogeneous case a = 1 for the case of arbitrary < a < 1. Our calculations 
predict that there occur a maximum of the relaxation rate constant at vary- 
ing this parameter. We carry out experimental measurements of the relaxation 
rate constant in porous glass Vycor and obtain data which agree well with our 
theoretical findings. 
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Table 1 



Characteristics of porous glasses. 
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Fig. 1. The dependence of the spectral density of the correlation function for spher- 
ical harmonics due to anomalous translational diffusion given by the (12) on the 
reduced Larmor frequency ULTa at different values of the parameter of inhomo- 
geneity a. The values of the latter parameter from the down line to the upper one 
respectively are: 1; 0.95; 0.9; 0.85. Here N is the number of spins in 1 cm^, d is the 
doubled radius of the molecules, t\ is the characteristic time for the homogeneous 
media given by the (16) and is the characteristic time for the inhomogeneous 
media given by the (15). 
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Fig. 2. The dependence of the contribution to the spin-lattice relaxation rate con- 
stant due to anomalous translational diffusion given by the (17) on the parameter 
of inhomogeneity a at different values of the reduced Larmor frequency ulTi. The 
values of the latter parameter from the down line to the upper one respectively are: 
1.; 0.4; 0.16; 0.06; 0.025; 0.01. Here N is the number of spins in 1 cm^, d is the 
doubled radius of the molecules, Ci is the diffusion coefficient for the homogeneous 
media given by the Stokes formula (14), ti is the characteristic time for the homo- 
geneous media given by the (16), 7 is the gyromagnetic ratio of the nucleus, I is 
their spin and h is the Planck constant. 
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Fig. 3. The experimental decay of the longitudinal magnetization for hexane in the 
porous glass Vycor 160 with the radius of the pores 160 A. Solid line is the result of 
the linear fit of the data. Here t is the current time of the magnetization's relaxation, 
Aq is the equilibrium value of the magnetization and A{t) is the current value of 
the magnetization. 
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Fig. 4. The experimental dependence of the spin-lattice relaxation rate constant 
on the reciprocal radius of the pores for the porous glasses Vycor at the Larmor 
frequency 19.5 MHz. Here Ti is the spin-lattice relaxation time and p is the radius 
of the pores measured in Engstroms A. The solid line is an empirical fitted one "to 
guide the eye" . 



25 



